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AERONAUTIC SYMBOLS 

1 . FUNDAMENTAL AND DERIVED UNITS 


' - v, , — 



Symbol 

. . V 

Metric 

English; , ' 

Unit 

Ahbrevia- 
, tioh 

umt 

Abbrevia- 

tion 

Leng:th 

1 

meter.: L 

second 

m 

s 

foot (or mile) 

second (or hpur).. ,.. 

ft (Or mi) 
see (or .hr) 

lb' ^ 

Force 

F -- ' - 

weight of 1 kilogram... 

weight of 1 pound. _i . . 

Power^i—^— 
Speed 

P 

. ■ V 

horsepower (metric)'. . . _ . 
/kilometers per hour... — 
\metera per second 

kph 
■ mps 

horsepower. 

miles per hour. . 

feet per second 

fps 




2. GENERAL SYMBOLS 


- WeigLt=mjr / , _ 

Staadard acceleration of gra'\Tity=9.80665 m/s* 
or 32.1740 ft/sec* 

Mass=^ 

Moment of inertia=mi*. (Indicate apds of 
radius of gyration i by proper subscript.) 
Coefficient of viscosity 


» ■ Kinematic viscosity ' _ 

p' Density, (mass "per unit volume) ^ ; ^ 

Standard density of dry air, 0.12497 kg-m"*-s* at 15° G" 
and 760 mm; or 0.002378 lb-ft~* sec* ~ 

Specific weight of “standard” air, 1.2255 kg/m* or . 
0.07651 Ib/cu ft , " ; ■ 


*. AERODYNAMIC SYMBOLS 




■"■Area- ■ v ' . . ' 

Area of wing 

Gap " 

Span 
Chord 

Aspect ratio, ^ 

True air s^eed ' - 

Dynainici)ressure,'^pF* : 

, Lift, absolute coefficient Oz,=-^ 

, Drag,"absolute coefficient, 

Profile drag, absolute coefficient 

Induced drag, absolute coefficient Co, 

Parasite drag, absolute coefficient 


Do 


Angle of setting of wings (relative to thrust line)^ 
Angle of stabilizer setting (relative to thrust 
' Ime) I • ^ ' : " ? 

Q Resultant moment 

O Resultant angular velocity - v 

B Reynolds nximber, p^ where Z is a linear dimen- 

sion (e.g., for an airfoil of 1.0 ft chord, 100 mph, 
standard pressure at 15° C, the corresponding 
- Reynolds number is 935,400; or for an airfoE_ 

of l.p m chord, 100 mps, the corresponding 
Rejmolds number is 6,865)000) 
a Angle of attack; - 

€ Angle of downwash 

<*o Angle of attack, infinite aspect ratio 

a< ' Angle of attack, induced . , , ; 

cca Angle of attack, absolute (measured from zero- 

lift position) .. 

y Fhght-path angle 


O 


C ■ Cross-wind force, absolute coefficient 
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ISOLATED AND CASCADE AIRFOILS WITH PRESCRIBED VELOCITY DISTRIBUTION 

By Arthur W. Goldstein and Meyer Jerison 


SUMMARY 

An exact solution of the 'problem of designing an airfoil with 
a prescribed velocity distribution on the suction surface in a 
given uniform flow of an incompressible perfect fluid is obtained 
by replacing the boundary of the airfoil by vortices. By this 
device, a method of solution is developed that is applicable both 
to isolated airfoils and, to airfoils in cascade. The conformal 
transformation of the designed airfoil into a circle can then be 
obtained and the velocity distribution at any angle of attack 
computed. Numerical illustrations of the method are given for 
the airfoil in cascade. 

INTRODUCTION 

The problem of increasing the output per stage in axial- 
flow compressors and turbines involves the use of high- 
solidity (closely spaced blades) stages of highly cambered 
blades. In addition, the velocity distribution must be care- 
fully selected as a function of arc length along the airfoil 
(blade section) boundary in order to avoid flow separation or 
excessively high local velocities. 

Several methods are available for obtaining an airfoil with 
a prescribed velocity distribution. The methods that lead to 
theoretically exact results are based on conformal-mapping 
theory. (See references 1 and 2.) In reference 3, Mutterperl 
extends the method of conformal mapping to solve the 
problem of computing a cascade of airfoils with prescribed 
velocity distribution but, for cascades with closely spaced or 
highly cambered airfoils, this procedure becomes very cum- 
bersome. Approximate solutions have been obtained by 
placing singularities such as vortices, sources, and sinks in a 
uniform stream. The shape of sections of airfoils in cascade 
can also be computed by distributing such singularities 
periodically throughout the region of the cascade, as described 
by Ackeret (reference 4). 

Because these vortex methods are not exact, a method 
with the vortices on the boundaries of the cascade airfoils 
was developed. This method gives a theoretically exact solu- 
tion without the computation difficulties encountered in 
conformal-mapping methods for highly cambered airfoils or 
closely spaced cascades. Furthermore, for the same accuracy 
in computing the airfoil shape, this vortex method requires 
the computation of fewer points than the method of conformal 
mapping because these points may be arbitrarily placed on 
the airfoil. The method may be applied to isolated airfoils 
and to airfoils in cascade. For the cascade, the inflow and 
discharge velocities and a velocity distribution on the surface 


of an airfoil are given and the shape of the airfoil is deter- 
mined. In some cases, the spacing of the blades is pre- 
assigned, which places a condition on the assumed velocity 
distribution. Once the airfoil shape has been evolved, the 
velocity distribution may be computed for any angle of 
attack by the method described in appendix A. The method 
of this paper was developed at the NACA Cleveland labora- 
tory during 1946. 

THEORY 

OUTLINE OF METHOD 

In reference 5, it is demonstrated that the two-dimensional 
potential flow about a body in a uniform stream can be 
represented by substituting for the body a sheet of vortices 
along its boundary. The vortex strength per arc length at 
any point is equal to the magnitude of the velocity at that 
point. A proof of this relation for the case of the cascade is 
given in appendix B. The problem of finding a shape with 
a prescribed velocity distribution when placed in a stream 
can then be stated: Given a vortex distribution, to find a 
contour which satisfies the condition that it will be a stream- 
line in the flow field induced by the uniform flow and the 
vortices distributed on the contour. 

The procedure of finding the shape begins by choosing an 
approximate shape and distributing the vortices on it. The 
stream function of the flow induced by the vortices and the 
uniform stream is computed at points on the boundary of 
the assumed shape. If this stream function is constant, the 
assumed shape is correct. Variations of the stream function 
are a measure of the deviation of the assumed shape from 
the correct one. These variations are used to distort the 
original shape into a new shape whose stream function is 
more nearly constant. The process is repeated until the 
variations become negligible. In the process of shape adjust- 
ment, the velocity is altered on the pressure surface. 

DERIVATION OF EQUATIONS FOR THE STREAM FUNCTION 

Isolated airfoil. — The. complex or reflected velocity 
w'{z) (which is the derivative of the complex potential func- 
tion w{z)) induced at the point z—-x-\-iy by a vortex of 
strength k located at Zo=Xo+iyo is 


(A summary of the principal symbols used in this report is 
given in appendix C.) 
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The complex velocity w'(z) induced by a uniform stream 
with complex velocity wj and a distribution of vortex 
strength per unit length y{zo) along a curve with coordinates 


IS 


\ , I 1 C y(^o) 

w{z)=w. 


dso 


( 1 ) 


where dso is the element of arc length along the curve. The 
complex potential at the point s is the integral of ut'iz) with 
respect to z, namely 

w{z)-=ZWu' + ^^jy{Zo) log (2-2») dSo ( 2 ) 

From reference 1 (notation modified), 

y{Zo) dso=w'{Zo) dzo=dw{Zo)=dip{Zo)-\-i d\p{z^ 

where 

q> velocity potential, R{w{z)] 

\j/ stream function, I[w{z)] 

Wlien equation (2) is applied to obtain the complex 
potential function at any point z in the flow field, the in- 
tegration must be carried out along the boundary of the 
body. Because this curve is a streamline, dif>=Q and, there- 
fore, equation (2) becomes 

vj{z)=2Wu' j log (z—Zo) d<p{z„) (2a) 


The imaginary part of equation (2a) is the stream function 
at the point 2 , 

4^=^ — xV,,+yV^—^j log ^{x—XoY+{y~ycYd<p{z„) (.3) 
where 

Vy y-component of uniform stream velocity V 
Vi ,r-component of uniform stream velocity V 

It is convenient to use the arc length along the airfoil as 
a parameter. If {x,y) is a point on the airfoil boundary, then 
s will denote the arc length there; similarly, will denote 
the arc Icngtli at (a„, y„). The vortex at ■% on the airfoil 


influences the stream function at the point s on the airfoil. 
The stream function induced at (x,y) by a vortex of unit 
strength at (xo, yo) is 


/i( 2 , Zo)=^ log [{x—XyY+ ly—VoY] (4) 

A plot in the {x,y) plane of curves for a constant /i ( 2 , 2 „) 
consists of concentric circles with center at {x„, y„). 

The velocity at the point So on the airfoil is the directional 
derivative <p'{so) of the potential along the streamline. 
If the velocity along the airfoil has been specified and an 
airfoil shape, has been assumed, the resultant stream func- 
tion along the boundary of the airfoil can be approximated 
by using the approximate shape in evaluating the integral 

'A(s) = '/'«(«) - ('b ■«») 'p' (s») dso (5) 

where 

tl'.Xs) stream f unction at (.r,?/) due to uniform stream, —xVyVyVi 
I total arc length of airfoil 

All variables are expressed in terms of the aredength para- 
meters s and s„. The integral in equation (5) can be evalu- 
ated either numerically or graphicallj" over the entire range 
of integration except in the region where a {^s—So) is small, 
for in this region /i(s,.s„) becomes infinite. This portion of 
the integral can be evaluated by approximating the airfoil 
boundary by a line segment. Then, 

So) log {s — So)- 


The prescribed velocity can be given in this region, which 
may be defined by s— a^So^s-fa, by a Taylor’s series as 
a function of about the point s. 


<p'{s^= (.S') -)- y" (■!>•) (■?„ 




where the primes indicate derivatives with respect to s. 
The integral is then 



^0) ^ (^*0) dSo~ 


rri [v.'(.so-by"(.v)(.-„- 

ij^a^o'(s) (log a — 1)+— (log 


*■)"!“•■■] ds'o 



(6) 


In most cases, only the first term need be used in equation (6). The same type of approximation can be used to evaluate 
a portion of the integral if the opposite side of the airfoil comes in the neighborhood of the point (/,?/). 

A more general equation applicable to a segment that does not pass through s is; 


1 rp+<^ . 


f 

'■Jjj 




-yo)V(.So) d(So)=^ \<p'{p) c log log {h^+h-)-2{c-h)V2h tan ‘ + 


’{'P) 


2 ! 


P K^'+c^) log {h^+c‘)-{k^+¥) log {h^+h’^)-{Y-h^)]V 


[c^ log (A^+O-6^ log (c®-6®) + 

2h\c^b)~2h? + - ■ • j (6a) 
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where h is the perpendicular distance from s to the segment, 
s„=p locates the foot of the segment, (p+6) and (p+c) 
are the limits of the integration of and approximately, 

<p"{v) = <p''{V^b)~h^'''{p+b) 

Equation (6a) may be used when the line segment is not 
of equal lengths on either side of the perpendicular foot or 
when <p'{s) or its derivatives are discontinuous at either 
(p + 6) or (p+c). If a=c= — b and h— 0 , equation (6a) 
reduces to equation (6). The size of a, 6, or c is deter- 
mined by the requirements that the segment in question 
be nearly straight (the approximation is of the second degree) 
and that <p'(so) be accurately represented by a Taylor’s 
series expansion of few terms. 



Ex it-velocify diagram Compietd-'y^iocity diagram 



Entrance -velocity diagram xoE Wm uii' 

Ref leaf ed- velocity diagram 

Figure 1. — Notation for cascade flow. 

Airfoils in cascade. — The expression for the complex po- 
tential for the flow about a cascade of airfoils is derived in 
appendix B. The notation is defined in figure 1. The 
equation that corresponds to equation (2a) for isolated air- 
foils is for a cascade of airfoils 

W(2) = 2W,„'-1-^ Jl0gj^sing(s — 2„)^ C?<p(2o) (7) 

where 

Wm mean stream velocity, which is one-half the sum of 
complex (reflected) velocities upstream and down- 
stream of cascade, — iVy 

S distance between successive airfoils in cascade 
The mean velocity Wm' corresponds to the uniform velocity 
Wy' of the isolated airfoil flow. 

The term zwj is the complex potential function resulting 
from the mean flow. In the integral, the element dip indi- 


cates the vortex-element strength and log [sin (tt/S) (z-Za)] 
represents the complex potential at the point 2 caused by an 
infinite row of unit vortices at Zo±nS where n=0, 1, 
2, . . .. The imaginary part of equation (7) is the stream 
function, 

J *So=l 

f2(.s,So) d<p{8o) ( 8 ) 

where 

f 2 (s,So) log j^sin^ ^ (x— a:„) -hsinh^ ^ (l/— 2/o) J 

is expressed in arc-length parameters and ^m(s) is the stream 
function at (x,y) induced by a mean stream whose complex 
velocity is Wm ' ; that is. 

The values of {x-Xo)IS and {y-yo)/S for various values of /2 
are given in table I. A plot of x-x„ and y-y„ for constant 
values 01 / 2 ( 2 , 2 ,) is shown in figure 2. These curves may be 


fi(x-Xy, y-yo) (vvo) 



(x-Xo) 


.2 


Figure 2.— Plot of curves for constant hix—xo, y—yo). 


interpreted as the streamlines of the flow induced by an 
infinite row of vortices of unit strength located at the points 
{xo±nS, yo), where n=Q, 1, 2, . . .. In the region of a 

vortex, the streamlines are nearly circles; that is, the flow is 
nearly that induced by an isolated vortex. At a distance 
from the vortex row, the streamlines are parallel lines, as in 
the flow pattern induced by a continuous uniform distribu- 
tion of vorticity along a straight line instead of a row of 
discrete vortices. The velocities on the two sides of such a 
vortex line are of equal magnitude but opposite in direction. 
This behavior of /2 for large \y — and also for small 

— — ^ can be described as follows; When both 

O ^ 

{x—Xo)jS and {y—y<>)IS are small, 

/2(2, z») log ;|2 [idi—^oY+(y-yoy] 




(9) 


which differs from /i (2, 2„) only by a constant. For large 
values of \y—yo\IS, irrespective of (x— x„)/iS and a constant 
term, 

/2(2,2,)«^3^ (10) 


v/hich is the stream function of a uniform stream parallel 
to the .T-axis. 
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Equation (8) can be used for computing the stream function along the boundary of an airfoil in cascade just as 
equation (5) is used for the isolated airfoil. The integral over the range in the neighborhood of the point s is obtained by 
using equation (9) for /2 (s,So). The result, derived in the sdme manner as equation (6), is 


V i2(s,s„)^'(s<,) ds<,= |a/(s) log a^— 1 ^ - g f- |^Iog ' ' ' | 


The more general equation (6a) is modified for cascades by 
multiplying the argument of all logarithms by the factor 
7t7*S-. 

ADJUSTMENT OF SHAPE 

If the stream function for the assumed airfoil has been 
computed and has been found to vary, the shape must then 
be adjusted to give a more nearly constant stream function. 
The shape changes are made by rotation of the body plus 
displacement of the individual points normal to the mean 
stream. The rotation is used to place the front stagnation 
point properly. 

Eotation of the airfoil. — In the formula for computing the 
stream function of an isolated airfoil, the contribution of a 
vortex element at {xg, yo) to the stream function of a point at 
(r, y) is dependent merely on the distance between the two 
points. Consequently, if the entire airfoil is rotated, the 
effect of the boundary vortices on the stream function at any 
point on the airfoil boundary will not change. The effect of 
the blade rotation on the stream function along the boundary 
is therefore determined by the change in relative position of 
the points in the uniform stream. The first adjustment in 
shape is a rigid rotation of the airfoil in order to obtain a 
more nearly constant stream function along its boundary. 

If the airfoil is rotated tlu’ough an angle /3, the stream 
function is so changed that (s) is a function of /3 and s and 


(6b) 

may be written ^ (s, /J). When /8=0, ^ (s, 0) is the original 
stream function before rotation. After rotation the new 
stream function \p (s, j8) may be expanded in a Taylor’s series 
about the point /3=0, 

^ 

Only the first two terms in this series will be used because (3 
is assumed to be small. The angle |8 is to be determined for 
tbe minimum mean-square deviation of the stream function 
from its mean value. Because the object of the rotation is 
essentially to adjust the shape of the nose, the rotation might 
also be made to reduce the root-mean-square deviation of the 
stream function to a minimum for a portion of the shape 
including the nose. 

The mean value of the stream function at any angle /3 is 

The difference between the new stream function tj/ (s, /3) and 
its mean value f (jS) is squared and integrated to obtain a 
measure of the variation of \p (s, /?) from the mean value at 
the new angle. The condition for obtaining a minimum 
root-mean-square deviation by adjusting ^ is 


-Jo' + * <‘2) 


The second integral vanishes by virtue of equation (11), 
which may also be used to eliminate yp (|3) from the remaining 
term. The solution for is 


«, 


- (*l 

yp{s) ds 
L-*./ 0 

[J 


1 

1 


^0 0/3 J Jo 

d^{sfi)~ 

- - 

ds 


In order to apply equation (14), d\p/d$ must be known at 
points along the boundary of the airfoil. For the isolated 
airfoil, the contribution of the vortices is unaffected by the 
rotation and therefore 

dyp dyp ji d , . y, dx . y- dy . 

If the airfoil is rotated about the point {Xc, y^), equation (15) 
becomes 


^=cos 13 [(x-Xc)V,+ (y-yg)V,] + 

sin /3 [{x—Xc) V,j— {y~yc) FJ (16) 

where (x, y) are the coordinates of the point before rotation. 
For small values of /3, equation (16) reduces to 

g=(.r-x,)F,-f(2/-2/,)F, (17) 

The choice of (x^, yc) will have no effect on the results in this 
case. 

When the airfoil in cascade is rotated, the change in the 
position of the vortices of the adjacent blade must be con- 
sidered. For the isolated airfoil, it was unnecessary to con- 
sider the change in position of the vortices because the 
influence of a vortex (equations (3), (4), and (5)) depended 
on the function/i, which is constant on circles. The influence 
of the vortices on the airfoil is therefore independent of 
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direction. Because the /2 contours are not circles, the rota- 
tion in cascade does have an effect, which is approximated 
by considering all closed /2 curves (/2<0) as circles in order 
that the effect of all vortices in the region ./2<C0 nisiy be 
neglected during rotation. The effect of all vortices in the 
region /2>0 is estimated by assuming that all the/2 contours 
for /2>0 are straight lines uniformly spaced. The flow 
coiTesponds to that between two infinite straight parallel 
vortex sheets of uniform strength per unit length. This 
flow induced by the vortices in the region /2>0 is in the 
a;-direction, and the direction of the flow induced by the 
vortices for which Vo^V is opposite in sense to that induced 
by the vortices for which yo<C.y- 

As the point being considered is changed, the regions for 
/2>0, yo^y, and /2^0, yo<^y will include different sections 
of the blades, and hence different vorticity, with the result 
that the a;-velocity component Vj: induced by the vortex 
sheets will vary with the point under consideration. The 
algebraic sum of the i^-component of the uniform flow 
velocity and the variable r-velocity I'j: induced by the vortices 
in the region /2>0 is to bo used like the velocity component 
in rotation of the isolated airfoil (equation (17)). The 
quantity in equation (17) is replaced by the' corresponding 
F,.,r= Fi-f The vortex strength per unit length at any 
point on the airfoil is equal to and, therefore, from 

equation (10) the r-coinponcnt of the velocity induced by the 

vortices is ‘p'(Sa) ds„, where the integration is carried 

out over the portion of the airfoil where Jo(s,s„)'^ 0 . A 
distinction must be made between the two regions yo<^y and 
1 /o^y because the induced velocity components have oppo- 
site dii'ections. 

The computed result of rotating an airfoil in cascade de 
pends upon the choice of (/cl/c)- In order to minimize the 
error involved, values of d 4 'l<W are reduced by choosing 
(rc,?/c) as the centroid of the vortex distribution on the airfoil. 
If the improvement in the mean-square deviation of ^ is 
small compared with its original value, it may be preferable 
to omit the rotation of the airfoil because of the error inherent 
in the approximation for dtpldfi. The decision should be made 
chiefly on how varies at the airfoil nose and whether it is 
approaching a constant value in this region with successive 
corrections of the shape. 

Distortion of the shape. — The stream function computed 
after the isolated airfoil has been rotated will, in general, 
still vary along the boundary. This variation can be re- 
duced by distorting the shape of the airfoil. If the distor- 
tion is small, the change in distance between any two points 
on the boundary will be small, although the change in the 
direction of a segment joining those points may be consider- 
able. The effect of the distortion on the contribution to 
the stream function of the vortices on the boundary is 
consequently neglected. The largest effect of the distortion 
will be to change the position of the boundary points in the 
uniform stream. The airfoil is therefore distorted in such a 
manner that the change in the contribution of the uniform 
stream to the stream function will eliminate the variations 
in stream function. For points directly opposite each other 
on the airfoil, the change in distance will be. of the same order 
of magnitude as the distortion. Consequently, distortions 


that result in change of thicloiess of the airfoil converge 
very slowly because of the inaccuracy of the fundamental 
assumption on which the correction is based. 

Thus, when the stream function along the boundary of 
the isolated airfoil is known, some number is arbitrarily 
chosen as the desired constant value of the stream function. 
If is the difference between the computed stream 

function at a point and the desu’ed constant, the point is 
moved a distance —Axl/fV perpendicular to the direction 
of the mean stream, where the direction of increasing uni- 
form stream function is taken as positive. The airfoil in a 
cascade is distorted in the same manner, by using the varying 
resultant local mean stream velocity -yjVx,/+V/; corrections 
are made with t/- equal to the mean value of \p on the airfoil. 
COMPUTATIONAL PROCEDURE FOR CASCADES 

CHOICE OF VELOCITY DISTRIBUTION 

Several factors influence the choice of the velocity distri- 
bution for which an airfoil is to be found. Especially in 
rotors, sturdy blades are required. Long thin tail sections 
must be avoided and where high rotative speeds and stresses 
occur, overhang of thin sections is likely to induce blade 
failure. The radial distribution of airfoil cross-sectional 
area is also fundamental in determining the blade-root 
stresses. Overhang can be reduced by proper choice of the 
velocity diagrams for the sections, but the other factors arc 
influenced chiefly by the thickness of the section. 

The desired thickness may be attained by first assuming a 
blade shape and spacing and by then using the stream-filament 
method of reference 6 to compute the velocity distribution 
over a portion of the airfoil that determines the thickness. 
The spacing may be regarded as fi.xed but the curvature can 
1)0 adjusted if local velocities arc too high for the desired 
thicl ■mess. This computed velocity will then serve as a 
guide to the choice of an airfoil velocity distribution, which 
sliould be chosen to avoid high velocity peaks and steep 
negative gradients. If the average of the velocities on 
opposite sides of the blade camber line is retained in the 
modification of the velocity distribution computed from the 
stream-filament method, the thickness will also be retained. 

Because of the irrotationality of the fluid motion, the 
velocity integral or circulation around the airfoil must be 
eciual to that around a blade but over a width equal to one 
blade space. Therefore, 

J" <p'{s) ds=r = iSXF„,i — Tx,2) 

where 

r circulation about airfoil 

Vx,i tangential velocity entering cascade 

tangential velocity leaving cascade 
This relation places a condition on the assumed velocity 
distribution. 

If the computations thus far have been made in order to 
select a velocity distribution for the airfoil cascade in a com- 
pressible fluid flow, an equivalent velocity distribution for 
the flow of an incompressible fluid must be determined 
before the blade shape can be computed by any method 
based on incompressible-flow tlieoiy. For subcritical flows, 
the directions of the inflow and discharge velocities are 
nearly the same for compressible and incompressible flows. 
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but for incompressible flow the component normal to cascade 
axis is the same upstream and downstream. The Karman- 
Tsien compressibility correction (reference 7) or that of 
Garrick and Kaplan (reference 8) may be applied to the 
velocity on the blade surface to estimate roughly the corres- 
ponding incompressible-flow velocity distribution. The re- 
sulting velocity distribution in any case must satisfy the 
circulation condition. This procedure does not give an 
exact solution for compressible flows, but the resultant com- 
pressible flow will have approximately the desired char- 
acteristics of low^ pressure gradients and no high velocity 
peaks. 

COMPUTATION OF AIKFOIL SHAPE FROM THE CHOSEN VELOCITY 
DISTRIBUTION 

The numerical computation of the quantities involved in 
the preceding analysis, particularly the function fo, is ex- 
tremely laborious when tables of/2(s,s„) are used. Most of 
the computations are therefore executed graphically. In the 
cascade example, the air was assumed to enter the cascade 
at an angle of 45° from the cascade axis and to leave at an 
angle of —30° from the cascade axis. The prescribed veloc- 
ity distribution is given in figure 3(a). The value of the lift 
coefficient for this airfoil is 3.1. The shapes of the isolated 
airfoil and the airfoil in cascade are computed by the fol- 
lowing steps: 



V 



(a) Initial airfoil. 

(b) Final airfoil. 


Fiol'RE 3.— rroscribed volocity distribution for thick airfoil in cascade. 

1. Curves for constant /i for the isolated airfoil, or con- 
stant /2 (fig. 2) for the airfoil in cascade, are drawn. This 
diagram should be made on some transparent material that 
will change neither in size nor shape. The coordinates of the 
curves for constant/2 ai'e given in table I. 

2. A desired velocity <p'(s) is chosen as a function of the 
arc length of the airfoil (fig. 3(a)). An airfoil shape having 
the desired total arc length is assumed and is drawn to the 
same scale as the plot of /i or fs. The drawing is made on 
grid paper and, in the case of a cascade, the x-axis coincides 
with the cascade axis (fig. 4). 

3. The velocity distribution (p'(s) is integi’ated to obtain 
the velocity potential <p(s). This function is plotted on the 


same chart as the assumed airfoil shape for the corresponding 
j/-coordinate, as shown in figure 4, by plotting both <p and 



Figl ke 4. — Plot of airfoil and velocitj' potential for use in computation. 


the 2/-coordinate of the airfoil against s on a supplementaiy 
graph. In regions of the airfoil where y varies little with s, 
that is, where the airfoil boundaiy is parallel to the x-direction, 
<p should be plotted against x in the same manner, as shown 
in figure 4. 

4. In order to find the stream function at a point (x,y) 
on the airfoil, /2 (s,s„) must be plotted as a function of ip(Sa) 
to evaluate the quantity J'Ji(s,So) dtp{s„) of equation (8). 
If the chart of /; is superimposed on the airfoil with one vor- 
tex center overlying the point {x,y), the value of/2 may be 
read at (x„,y„) and the corresponding value of <p{x„,yo) may 
also be read from the plot of tp{Xo,yo)- The value of /2(.s,s«) 
is the same as would have been obtained by centering the 
chart on (r„,y„) because of the symmetry of the function. 
A succession of values of 9? and/2 are obtained in this fashion 
for various positions {xo,yo) that intersect the contours, 
and a plot of these points (/2,<p) may be made for the assumed 
position (x,y). This procedure is illustrated in figure 5 for 
a particular point {x,y) on which the chart is centered. 
The readings for a particular (xo,?/„) are shown by the arrowed 
lines. The points 1 to 6 on the blade are shown on the 
corresponding /2 curve. The discontinuity of between 
points 1 and 6 is the circulation. The discontinuity between 

4 and 5 represents the region where ,/ approaches — «> . 

5. The proper method of integration then proceeds from 
1 through 6 to 7 and then to the origin, with constant /:> 
from 4 to 5. The region from 4 to 5 witli /o approaching 
— CD is computed by equation (6) or (fib); the constant a is 
assumed to bo the radius of the near-circle, which corres- 
ponds to tlie value of /2 where the discontinuity from 4 to 

5 occurs. 

The total area including this small addition is 

J ‘p'(.So)Ms,So) ds= j Ms,s„) d<p 

which is the stream function due to vortices on the entire set 
of airfoils in cascade. Wliere /2=0 at the points A, B, C, 
and D (fig. 5), the values of <p are noted as <pa(s), <pb(s), ‘Pc(s), 
and <rD(s). These values are used in computing the stream- 
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function change caused by rotating the blade. The stream 
function at the point {x, y) may now be computed from 
•equation (8) or (5), and 

^,n= — VyX+V^y 



Figure 5 —Superposition of figures 2 and 4 to obtain plot of /a against tp. 


A plot of the stream function (variation from the mean 
value) is shown in figure 6 for the initially assumed shape. 
Corresponding points on adjacent airfoils have a difference 
of Ai/'/'FyiS' equal to 1.0. 



Figure 6. — Variation in stream function along initial shape and first and seventh approxi- 
mations of airfoil cascade. 


6. When ^(s) is known at a sufficient number of points, 
the airfoil may be rotated as previously described. For the 
isolated airfoil, equations (14) and (17) may be used directly. 
For the airfoil in cascade, the coordinates of the centroid of 
the airfoil must first be computed by 



Before equation (17) can be used to compute d^/d/3, the vari- 
able quantity Vx.r must be computed. The vortices in the 
region/2)>0 are considered to be uniformly distributed along 
the cascade axis and the velocity induced by such a distri- 
bution is 


where y is the vortex strength per unit length along the cas- 
cade axis for/2>0. Therefore, 



where the integral is to be taken over the regions /2>0. The 
region /2)>0, 2/o>2/ contributed a positive component to Vx, 
whereas the region /2)>0, yo<^y contributes a negative com- 
ponent. The computation is simply carried out by making 
use of the fact that the integral for Vx is the difference between 
values of (p at points where /2=0. The values of ¥’a(«o), 
¥’b(so), ^c(So), and <pd{So) from step 5 are used at this point 
to obtain 


2vxS—J'<p'(So) dSo=(pA — ‘PD-^^—i<Pc—<PB) (18) 


where r is introduced because of the discontinuity in tp at 
the trailing edge. The sum cp^ — v’o+r gives the effect of 
the vorticity in the region ^(s, s»)]>0 near the trailing edge, 
and the term tpc — <pb gives the effect of the vorticity in the 
region/2 (s,s„)>0 near the leading edge. If either the leading 
edge or the trailing edge lies in the region /2(s,s„)<(0, only 
two points of intersection will remain and one of the two 
groups of terms in equation (18) will vanish. The quantity 

f dSo is added to the j;-component of the original 

2jo j 


uniform stream velocity and the quantity d^l/jd^ of equa- 
tion (17) may be computed for a number of points and the 
angle j8 computed from equation (14), using the values of {Xc,yc) 
just determined. After these computations have been made, 
the airfoil is rotated through the angle |8, and the value 


is assigned as the value of the stream function of 
dp 

the point after rotation. 

7. A value of ^(s) is Imown at points along the airfoil 
boundary. The mean value over the airfoil is subtracted 
from yj/ leaving A^. For the isolated airfoil, no subtraction 


is necessary. 


Each point is moved a distance 


M/ 

^jVx/+Vy^ 


in the direction perpendicular to the velocity computed in 
step 6. The curve joining the points in their new positions 
is the adjusted airfoil. 

8. The total arc length of the adjusted airfoil will be 
different from the original one, in general, although local 
changes in length will be negligible. The airfoil is so scaled 
that the length of the suction side is the same length as it was 
before distortion because this surface is the critical surface 
of the airfoil. This process will result in a change in length 
of the pressure side. The velocity over the pressure side <p' (s) 


821313- 
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must then be altered in such a manner that the difference 
in potential between the two stagnation points remains the 
same. As a result, the quantities that retain specified values 
are the length and the velocity distribution on the suction 
side and the circulation around the airfoil. The entire 
procedure is repeated with the adjusted shape until the 
variations in the stream function result in very little change 
in the shape of the airfoil. 

DISCUSSION OF EXAMPLES AND TECHNIQUES 

For the example being computed, the stream functions 
obtained for the initially assumed shape and the first and 
seventh approximations are plotted against the arc length 
(fig. 6), which is taken as zero at the trailing edge and pro- 
ceeds counterclockwise around the airfoil as shown in 
figure 7. I'he fact that for the initial shape is positive 
over the first half of the arc length and negative over the 
second half indicates that it is too thick because the required 
distortion in shape, will make it thinner. The change in 
thickness results in a change in velocity distribution over the 
pressure side of the airfoil in order to maintain the desired 
circulation. The velocity that was originally assumed, 
which is equal to the vorticity per unit length distributed on 

l^j -f-jojl Qfp fQ // 

Final airfaH 



Figure 7.— Initial shape and final approximation of thick airfoil showing cascade spacing. 

the initial airfoil, is shown in figure 3 (a) and the velocity 
over the final shape in figure 3 (b). The length of the 
pressure side has increased and the velocity has decreased 
in the portion of 1:1.1. 


Over the section of the aufoil that has collapsed to zero 
thickness, the surface velocities of figure 3 (b) may not have 
been obtained, but the loading (circulation per unit arc 
length), which is the difference in the velocities on opposite 
sides, has been realized. In practice, this coUapse is pre- 
vented by increasing the assumed velocity on the airfoil 
surface. 

If the initially assumed airfoil shape has a thickness that 
differs considerably from the correct one, the process of shape 
adjustment will converge rather slowly. The labor can be 
reduced, however, by computing the stream function at a 
few points on the airfoil and locating these points to deter- 
mine the thickness. This procedure is followed for the first 
few approximations until the thiclmess of the airfoil is fairly 
accurate. The stream function is then computed at a 
larger number of points, particularly near the leading edge, 
in order to get more detail of the shape. 

Arbitrary specification of a velocity distribution may 
result, not in a physically real airfoil, but in a figure-8 
shape or a collapsed shape (zero thickness over a portion of 
the blade) . The velocity distribution must then be modified 
to obtain a real shape ; these modifications should be selected 
to keep the desirable properties of the original distribution. 
Velocity peaks and steep velocity gradients, which tend to 
occur on the suction side of an airfoil, are to be avoided. 
If the airfoil collapses, the vorticities of the two sides tend to 
cancel each other and the remaining vorticity represents the 
difference in velocity across the thin airfoil rather than the 
velocity along the boundary. 

The method was also applied to the design of a thin airfoil 
(camber line) in a cascade. The vortex distribution is 
equivalent to load distribution (difference in velocity across 
the airfoil) rather than velocity as in the case of a thick 
airfoil. The velocity diagram for the cascade and the desired 
load distribution for the thin airfoil are shown in figure 8. 
The value of the lift coefficient of the resultant airfoil is 4.1. 



s 

.S’ 

Figure 8. — ^Velocity diagram for cascade and prescribed load distribution for thin airfoil in 

cascade. 
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isolated and cascade aieeoils with prescribed velocity distribution 


The initial shape was obtained by assuming zero spacing 
between the airfoils. The initial shape and the first and 
third approximations to the airfoil shape are shown in figure 9. 

Initial airfoil 

First approximation 

Third approxima'i'ion 




Figure 9. — Assumed shape and first and third approximations of thin airfoil showing cascade 

spacing. 


The second and third approximations differ very little. 
The third approximation is redrawn in this diagram to show 


the spacing between airfoils. The convergence of the method 
is shown graphically m figure 10. The variation A^' of the 



Fiqube 10. — Variation in stream function for successive approximation of thin airfoil in 

cascade. 


stream function from its mean is divided by VS to make it 
dimensionless and is plotted against the arc length along the 
airfoil where s=0 at the trailing edge. The stream function 
computed on the second approximation is nearly constant, 
which gives the third approximation almost the same shape 
as the second one. The rapid adjustment of camber con- 
trasts with the slow adjustrnent of thickness. 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 4, 1947. 




APPENDIX A 

VELOCITY DISTRIBUTION ON THE DERIVED AIRFOIL AT DIFFERENT FLOW ANGLES 


Conformal mapping. — When an airfoil is given, the ve- 
locity distribution over its surface must frequently be found 
at different angles of attack. This problem may be solved 
by the method of conformal mapping, which consists in 
mapping the region exterior to the airfoil on the exterior of 
a circle. The velocity around the airfoil is obtained from the 
known velocity around the circle. Procedures for finding 
the function that maps a given airfoil into a circle are pre- 
sented in references 1 and 9 for the isolated airfoil and refer- 
ences 3 and 10 for the airfoil in cascade. 

In general, the procedure for finding the mapping function 
of an airfoil is a laborious one. But when, as in the present 
case, the velocity distribution over the airfoil at a particular 
angle of attack is known, the correspondence between points 
on the airfoil and on the circle, and hence the flow velocity 
at other angles of attack, can be obtained very easily. 
Indeed, the correspondence of points and the velocities for 
various angles of attack can be obtained by the method given 
in reference 11 from the initial data without knowing the 
airfoil shape, because the complex potentials of the airfoil 
plane and the mapping-circle plane are equal. Before the 
airfoil is designed it is therefore possible to check whether 
the airfoil to be computed will be satisfactorj’’ under condi- 
tions different from the design condition. 

Isolated airfoil. — The flow about any airfoil shape can be 
mapped on the flow about a unit circle in such a way that 
corresponding points have the same potential. The flow 
about the airfoil is given and the potential function <p(s) at 
each point is computed. If the potential function on the 
airfoil is computed by integrating the velocity from the 
stagnation point at the trailing edge in a counterclockwise 
direction around the airfoil oriented as in figure 1, the poten- 
tial will be zero at the trailing edge, decrease to a minimum 
<Pmin at the stagnation point at the leading edge, and then 
increase to a value equal to the circulation r at the trailing 
edge. The corresponding flow about the circle is determined 
by the conditions that it must have the same values of <pmtn 
and r for a correspondence to exist between all airfoil and 
circle points. If dr is the central angle of the stagnation 
point on the circle that corresponds to the trailing edge of 
the airfoil, 

^^^=-(cot0r+0r+x/2) (Al) 


Equation (Al) can be solved numerically for dx because all 
the other quantities are known. The velocity at infinity in 
the circle plane Vc can then be determined from the Kutta- 
Joukowsky condition, which requires that dx be a stagnation 
point; that is. 


F.= - 


r 

Itt sin dx 


(A2) 


The velocity potential at points on the circle is 


<Pc=—2Vc cos + cosdx—-^9x 


2tt "" 


(A3) 


The quantity 2Vc cos dx — ^®ris a constant that is subtracted 


in order to make at the stagnation point correspond- 
ing to the trailing edge. 

The correspondence of points on the airfoil with points on 
the circle is obtained by associating points where <p(s)=tpc. 
The velocity on the circle at a uniform stream flow angle a is 


v^{d,a)=2 Vc [sin (0 -f- a) — sin (e r a) ] ( A4) 


The nature of the conformal transformation is such that the 
ratio of the velocity at a point on the airfoil to the velocity 
at the corresponding point on the circle is independent of 
angle of attack. Therefore, the velocity (s) on the airfoil 
at flow angle a is 


Vc(0,a) Vc(e,0) 


(A5) 


where the design flow angle is taken as zero. Equation (A5) 
can be used to compute the velocity distribution on the air- 
foil except at the two points that were stagnation points at 
the design angle of attack. 

Airfoils in cascade. — The flow about a cascade of airfoils 
can be mapped conformally into the flow about a unit circle 
with two singular points located on the real axis symmetri- 
cally with respect to the center of the circle. These singular 
points correspond to the points at infinity in front of and 
behind the cascade, respectively. In a cascade of airfoils, 
the distance of these points from the center of the circle is 
uniquely determined by the same conditions that determine 
the flow about the circle in the isolated case; namely, the 
circulation per airfoil, the velocity potential at the leading 
edge, the blade spacing, and the upstream and downstream 
flow angles. 

The distance from the singular points to the center of the 
circle is denoted by e^. The flow about the circle is such 
that the location of the stagnation points ds is determined 
by the relation 


r sin 6s 
2VS~sm]i 


cos ds 
cosh K 


sin X 


(A6) 


where X is the angle of inclination of the mean stream to the 
normal to the cascade axis. (See reference 6 for details.) 
The quantities F, V, S, and X are known from the flow in 
the cascade plane and therefore equation (A6) provides a 
relation between K and the location of the stagnation points. 
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The velocity potential at anj'- point on the circle is 


<Pc,c=- 


vs r . 


(sin X .an- X ianh- )+L .an- 

-cos X tanh-' t»n-' 


rFfS / • , , _i sin Or 
(smXtan 


(A7) 


where 6t is the particular value of 0^ corresponding to the trailing edge of the airfoil. The expression in brackets is a 
constant so chosen that the potential will vanish at the stagnation point corresponding to the trailing edge of the airfoil. In 
order to map the cascade on the circle, K must be found so that the value of ipc.c at d^, the stagnation point corresponding 
to the leading edge of the airfoil, is equal to tpmtn, the value of the velocity potential there. The identity 


/ sin Bs Y 
Vsinh K) 


sinh^ 



cosh- K^l 


is used to eliminate 


cos Bs 
cosh K 


from equation (A6) to give 


sin Bs 
sinh K 


r 

2VS 


cosh^ K sin X ± cos 


X-y cosh^ K—cos^ ^~(^ 2Vs } ^ sinh^ K 

cosh^ AT— cos^ X 


(AS) 


In successive approximations, a value of K is assumed and equations (AS) and (A6) are used to find 


sin Bk cos 0JV 


sin Bt 


and 


cos B-r 


These values are inserted into equation (A7) to find (pc.c at B—Bn- If <pc,c{Bn) is not equal to (p„ 


Snh K’ cosh if’ vmues arc uiseibeu iiiio equation K.AI ) to nnu ipc,c 

another value of K is chosen, on the premise that <Pc.c{Bfr) will decrease as K is decreased. When iPc,c{6n) is evaluated, 
care should be taken to use consistent values of the inverse tangents. After two values of K and v’c.c(^iv) are determined, 
interpolation or extrapolation may be used for new values of K. 

A'V hen K has been found, it is used in equation (A7) to evaluate ^’c.c at values of B all around the circle. A point on 
the circle corresponds to the point on the airfoil where <p{s)=ipc,c- The velocity at the point B on the circle is 


VS sinh 2K 

X cosh 2i?— COS 20 



cos B 
cosh K 


cos Bt \ 
cosh KJ 


-f-sin X 


( sin B sin Bt 

sinh K sinh Kj_ 


(A9) 


and the velocity on the airfoil at any other mean 

flow angle X + a is 

^.'(i.')=*^c.c(0,X+a) (AlO) 

as in the case of the isolated airfoil. 

The designed airfoil was mapped on the unit circle by the 
method described. The constant K, the natural logarithm 
of the distance from the singular points to the center of the 
unit circle, is 0.075. The correspondence between points 
on the airfoil and those on the circle is plotted in figure 11, 
which shows the arc length of the airfoil as a function of the 
central angle of the circle. The velocity at any point on the 
airfoil for any angle of attack a may be obtained from 
equations (A9) and (AlO), the velocity distribution as in 
figure 3 (b), and the relation between s and B as in figui'e 11 . 

The ratio ^ is equal to dBlds (radians) and need be 

computed only once for any given airfoil. 



Figure 11.— Correspondence between points on airfoU and points on unit circle by conformal 

transformation. 


APPENDIX B 

DERIVATION OF THE CASCADE EQUATION 



An equation is to be developed for the complex velocity 
at any point in the field of flow of a fluid past a row of equally 
spaced, congruent bodies. Coordinates axes are chosen with 
the origin inside one of the bodies and the x-axis in the 
direction of the row. (See fig. 12.) The hody containing 



Figure 12.— Diagram for derivation of equation for flow about ca.seaclc. 


the origin is denoted by Bq, bodies along the positive di- 
rection of the a;-axis hy Bi, B 2 , . . ., and along the negative 
direction of the .r-axis by B-i, B^ 2 , .... A circle A of small 
radius is drawn about the point 2 where the velocity is to 
be determined. A rectangle R is drawn with its center at 
the origin and its sides parallel to the axes of length (2iV-l-l)jS' 
and width 2t, which contains the bodies B_n, . . . /?-i, 
Bo, Bi, . . . Btf, and the circle A. If a side of the rectangle 
intersects one of the bodies, the side may be distorted to go 
around the body with no essential change in the proof. 
The function w'(zo)l2o — s) is an analytical function of z in 
the region inside the rectangle R but outside the bodies 
and the circle A. 

Therefore 





dzo=0 (Bl) 


The first integral can be broken iqt into four integrals, 
one along each side of the rectangle, namely. 


r w'(x-it) n w’l(N+l l 2)S+iy„] , C-> 

Jr2„ — Z ° J -(ZV+1/2)S *0 — — 2 ° J_i (A^-|- 1/2) 2/0 — 2 ° J (W+1/2)S XoT it 2 ° ' Jt 


w'[-{N+ll2)S+iy o] . , 
~{N+l/2)S+iyo~2 


(B2) 


In an evaluation of these integrals, the function w'(zo) is 
periodic, with period S, and approaches a constant value 
infinitely far from the cascade; that is, 


and 


w' (x„+iyo)^W 2 ' as 
w'{x„+iyo)-^Wi' as 


From the last of these conditions, it follows that 


where 


w'{x„— it) = W 3 ' {Xo — it) + w, ’ 


Wz {xo-it)-^^ as t- 


Therefore, the first integral on the right side of equation 
(B2) is 


p(w+i/ 2 )s w'{x„—it) 
J-(iV+l/ 2 )S Xo — it — z 


pW+l,-2)S ^ pA’+l/2)S Wg'/.gp— if ) 

J -(X+H 2 )SXo—it — 2 J-(X+U 2 )S Xo — it — z 


(B3) 


The first of these integrals is 




(iv+i/2)s [(;V'’+l/2)-S'-if— z] . , 

(i\'+l/ 2 )S — d—z [—{A-{-l/2)o — it — 4.J 


asAi-»co andf-^=o, provided that f/(A^»S/) ^0. The last integral in equation (B3) is 


r 


w/(Xo-it ) 
(N+1J2)S Xo — it — 2 


dx 


(»+>f2)s Wz'{Xo-it) 

-1/2)S Xo — it — Z 


„= z: f 

n=~N J (n 

_ z wo'jxo-it) , 

n=-x J-s /2 XoAnS — it — 2 ° 

J-s/2 Xo — It — 2 n=l J-S/2 


2{Xo—it — z)wz {x,Q—it) dxp 
{Xo—it—z)^— rdS^ 


i 


4 . 
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If t is chosen sufficiently large so that| W 3 '(a;o— ii) | ■<«, where e is any preassigned positive number, the moduli of the integrals 
are less than 


r dx, N pS/2 2\x„-it-z\dXo 1 < r 

U-s /2 \Xo—it—zy hi J-s /2 \{x„—it-zY—n^S^\\ -* LJ- 


dxo 

sn -yJiXo—xY+it+y) 


-u'h 2 V {Xq — x)^ H~ dxi> ~j 

’■ hi J-s /2 (x„—xY+(t+yY—'>^^^^-i 


When iV-^ , this quantity approaches ^ ^ 

rs /2 _ r„ , ~1 conditions, 

« -sCot (v+i/as 'i- -a 

J(Ar+i/ 2 )s Xo-j-it—s 

This integral is finite and, because e can be made arbitraxdly The second and fourth integrals on the right side of equa- 

small as <» ^ the last integral in equation (B3) approaches tion (B2) can be evaluated by combining them. Because 


* / s /2 S \3^ (a;^— ar)H- (<+2/)^J dx^ 


zero. Therefore, 


p(N+l/2)S w'{Xo-it) , 
j-(N+mS Xa-it-Z 


w' is periodic, 

ia'[(iV+ \l2)S+iy,]=w'[- (iV+ l/2)S+iyJ 
and therefore. 


n w' [(N+l/2)S+iy„] . r-‘ w'[-(N+l/2)S+iy,] _ -2(N+i/2)Sw'l(N+l/2)S+iy.] . , 

(iV+l/2)S+'i2/,-2"'^^“+J, -(N+l/2)S+iy,-z^^^‘' J-t (iy,-zy—(^+l/2YS^ 

The velocity w'[(iV+ l/2)»S'+ii/„] is bounded for all values of 2 /„; that is, there is a constant IF such that 
|w'[(A7'+l/2)S+i2/o]Kl'T- The absolute value of the integral is less than 


2S(iV+l/2)IF J* 


, \{iy„-zY-{N^ll2YS^\ 


;2S(.v+,/2)ir£^ 


dy^ 

-2/)^+(iV+l/2)2S^-x2 


2^(W+1/2)IF -t—y 

'V(W+l/2")^5'^-r^L v(iV+l/2)"5'^-a:^ v(^+ 


As CO and this quantity approaches zero. It has ^y the residue theorem, 


been shown, therefore, that when and >0, 


r w'{z„) 

}aZo — Z 


dz„—2Tviw' {z) 


f d2„->vi (Wj' + w/) ( 

JR *0 — 2 

± r d2,= i: , f 

n — —NjBn ^ u=:—Y\jB'i'^o i 


(B4) periodicity of w' {z) implies that 


= r %ihdz„+iz\ 


°^Xs |w'(2„) cot I {Zo—Z) 


as CO . I The complex potential is obtained from equation (B7) by 

When equations (B4), (B5), and (B6) are substituted into integrating with respect to £ and neglecting the arbitrary 

equation (Bl), the expression for the complex velocity is constant, 

obtained: y,( 2 ) = zwj + ^ w'{Zo) log sin ^ ( 2 — 2 „) dzo (B8) 

{z)-^{wY + wY) 2 ^. L ^w'{z„) cot ^ (Zo—z) dZo (B7) I where Wm' = ^' is the mean stream velocity. 



APPENDIX C 

SYMBOLS 


The principal symbols used thi’oughout the report are 
listed here for convenience of reference. 

/i = 4 ^ log {^—^oY+iv—yo ) -] 

J 2 = 4 “ log sin^ ^ (a;-a;„) + sinh- ^ (?/— 2 /„)J 

K natural logarithm of distance from singular point to 
center of circle corresponding to cascade airfoil 
I total arc length of airfoil 

8 distance between successive airfoils in cascade 
s arc-length parameter corresponding to z 

So arc-length parameter corresponding to z„ 

V magnitude of uniform or mean stream velocity in air- 
foil or cascade plane (fig. 1) 

Vo magnitude of uniform stream velocity in circle plane 
Vx a:-component of uniform or mean stream velocity F 
r resultant local mean stream a;-component of velocity V 
Vy y-component of uniform or mean stream velocity F 
Vo local velocity on circle corresponding to isolated 
airfoil 

Vo,c local velocity on circle cori’esponding to airfoil in 
cascade 

velocity induced by vortices in region / 2>0 
w complex potential function, <p+i\f/ 
w„' complex velocity of mean stream for airfoil in cascade 

{Wl' + W2) = Vx — iV^ 

wj complex velocity of uniform stream for isolated airfoil, 

b X' y 

X real part of z 

^c.Vc coordinates of point about which airfoil is rotated 
(centroid of vortex distribution for cascade airfoils) 
y imaginary part of z 

s coordinate of point where stream function is com- 
puted, x-{-iy 

Zo coordinate of point where vortex is located, Xo+iyo 
a angle of inclination of uniform stream velocity to 
x-axis 

/3 angle through which airfoil is rotated 

r circulation about airfoil 

-/(Zo) vortex strength per unit arc length at z„ 

0 central angle of circle 

angle of stagnation point on circle corresponding to 
leading edge of airfoil 

St angle of stagnation point on circle corresponding to 
trailing edge of airfoil 


X angle of inclination of mean flow to normal to cascade 
axis (fig. 1) 

<p velocity potential on airfoil, B[w{z)] 

<Pa,<Pb, values of tp at points A, B, C, D, respectively, where 
>Pc><pD curve of ^(s„) intersects /a (s, So)=0 (See fig. 5.) 

(Pc velocity potential on circle corresponding to isolated 
airfoil 

^c.c velocity potential on circle corresponding to airfoil in 
cascade 

<Pmin velocity potential at leading edge of airfoil 

^ stream function, I[w(z)] 

\!/m stream function of mean stream of cascade flow 

stream function of uniform stream flowing about 
isolated .airfoil 

\l/ mean value of stream function over airfoil 

Axp variation of stream function, i/'— i/' 

Subscripts 1 and 2 when appended to w' , V, and b’^. indicate 
inflow and discharge values, respectively. 
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ISOLATED AND CASCADE AIRFOILS WITH PRESCRIBED VELOCITY DISTRIBUTION 

TABLE I. COORDINATES OF /2(x-x„, y-yo) 

(a) Values of {y—y„)IS 


\ X—Xo 

\-^ 
h\ \ 

i \ 

0 

0. 05 

0. 10 

0.15 

0.20 

0. 25 

0. 30 

0. 35 

0. 40 

0.45 

0. 50 

-0. 40 

0. 0257 











-.38 

.0292 











-.36 

.0331 











-. 34 

.0376 











-.32 

' .0425 • 











-0. 30 

0. 0481 











-.28 

.0545 

0.0229 










-.26 

.0618 

.0371 










-.24 

.0699 

.0497 










-.22 

.0791 

.0621 










-0. 20 

0. 0894 

0. 0750 










-.18 

.1010 

.0887 

0. 0296 









-.16 

. 1140 

. 1035 

.0620 









-.14 

.1286 

.1195 

.0871 









-.12 

. 1447 

.1369 

.1107 

0.0392 








-0. 10 

0. 1626 

0. 1560 

0. 1344 

0. 0881 








-.08 

. 1824 

.1768 

.1588 

.1241 

0.0454 







-.06 

.2041 

.1993 

.1844 

.1572 

.1103 







-.04 

.2277 

.2236 

.2113 

.1896 

.1558 

0. 1014 






-.02 

.2532 

.2498 

.2396 

.2222 

.1066 

.1608 

0. 1096 





0 

0. 2805 

0. 2778 

0. 2094 

0. 2553 

0. 2354 

0. 2006 

0. 1777 

0. 1400 

0. 0969 

0. 0496 

0. 0000 

.02 

.3097 

.3074 

.3006 

.2892 

.2737 

.2542 

. 2318 

.2081 

.1858 

. 1692 

. 1629 

.04 

. 3405 

.3386 

.3331 

.3239 

.3117 

.2968 

. 2804 

. 2638 

. 2491 

. 2389 

.2352 

.06 

.3728 

.3713 

. 3668 

.3595 

.3498 

.3384 

.3260 

.3139 

.3036 

.2965 

.2940 

.08 

.4004 

.4052 

.4016 

. 3958 

.3881 

.3793 

.3698 

.3608 

.3533 

.3482 

.3464 

0. 10 

0.4412 

0.4402 

0. 4373 

0. 4327 

0. 4267 

0.4198 

0. 4126 

0.4058 

0.4001 

0. 3964 

0. 3951 

. 12 

. 4709 

. 4761 

. 4739 

.4702 

.4655 

. 4601 

. 4545 

. 4493 

.4451 

.4423 

.4413 

. 14 

.5135 

.5129 

.5111 

.5082 

.5045 

.5003 

. 4960 

. 4920 

.4888 

.4867 

. 4859 

. 16 

.5508 

.5503 

..54 88 

.5466 

.5437 

.5404 

.5371 

.5340 

. 5316 

.5299 

.5294 

. 18 

.5886 

. 5882 

.5871 

. 5853 

. 5830 

. 5805 

.5779 

. .5735 

. 5730 

. 5724 

. 5720 

0. 20 

0. 6269 

0, 6266 

0. 6257 

0. 6243 

0. 6225 

0. 6205 

0.6185 

0. 6167 

0. 6152 

0.6143 

0. « H 0 

.22 

. 6655 

. 0653 

. 61)46 

. 6635 

.6621 

. 0606 

.6590 

. 0576 

. 6565 

.6557 

. 6555 

. 24 

.7044 

.7042 

.7037 

.7029 

.7018 

. 7006 

.6994 

. 6983 

. 6974 

. 6968 

. 6966 

.26 

.7436 

.7434 

.7430 

.7424 

.7415 

. 7406 

.7396 

. 7588 

. 7381 

. 7377 

.7375 

.28 

.7829 

.7828 

. 7825 

.7820 

.7813 

.7806 

.7798 

.7792 

.7787 

.7783 

.7782 

0. 30 

0. 8224 

0. 8223 

0. 8221 

0. 8217 

0. 8211 

0. 8206 

0. 8200 

0.8195 

0.8191 

0.8188 

0.8187 

.32 

.8620 

.8619 

. 8617 

.8614 

.8610 

.8606 

.8601 

. 8597 

. 8594 

.8592 

.8592 

.34 

.9017 

.9016 

.9015 

.9012 

.9009 

.9006 

.9002 

. 8999 

. 8997 

.8995 

.8995 

.36 

.9415 

.9414 

.9413 

.9411 

.9409 

.9406 

.9403 

.9401 

. 9399 

.9398 

.9397 

.38 

.9813 

. 9812 

.9811 

.9810 

.9808 

. 9806 

. 9804 

. 9802 

. 9800 

.9800 

.9799 

0. 40 

1.0211 

1.0211 

1. 0210 

1.0209 

1.0207 

1. 0206 

1.0204 

1.0203 

1.0201 

1.0201 

1.0201 

.42 

1.0610 

1.0610 

1.0609 

1.0608 

1.0607 

1.0606 

1.0604 

1.0603 

1.0602 

1.0602 

1.0602 

.44 

1. 1009 

1. 1009 

1. 1008 

1. 1008 

1. 1007 

1. 1006 

1. 1005 

1. 1004 

1. 1003 

1. 1003 

1. 1003 

. 40 

1. 1408 

1. 1408 

1. 1408 

1. 1407 

1.1407 

1. 1406 

1. 1405 

1. 1404 

1. 1404 

1. 1403 

1. 1403 

.48 

1. 1808 

1. 1807 

1. 1807 

1.1807 

1. 1806 

1.1806 

1. 1805 

1. 1805 

1. 1804 

I . 1804 

1, 1804 

0. 50 

1. 2207 

1. 2207 

1. 2207 

1.2206 

1. 2206 

1. 2206 

1.2205 

1. 2205 

I . 2204 

1 2204 

1. 2204 

.52 

1. 2607 

1. 2607 

1. 2607 

1.2606 

1.2606 

1.2606 

1.2605 

1. 260.5 

1. 2605 

1 . 2005 

1 . 2605 

.54 

1.3007 

1.3006 

1.3006 

1.3006 

1.3006 

1 . 3006 

1.3005 

1.3005 

1. 3005 

1. 3005 

1.3005 

. 56 

1.3406 

1.3406 

1. 3406 

1.3406 

1.34( H ) 

1.3406 

1.3405 

1 . 3405 

1.3405 

1.3405 

1.3405 

. 58 

1.3806 

1. 3806 

1.3806 

3 . 3806 

1.3806 

1.3806 

1.3805 

1 . 3805 

1.3805 

1 . 3805 

1.3805 

0.60 

1.4206 

1.4206 

1.4206 

1.4206 

1.4206 

1. 4206 

1.4205 

1 . 4205 

1.4205 

1.4205 

1.4205 

.62 

1. 4606 

1.4606 

1.4606 

1.4006 

1.4606 

1.4606 

1.4605 

J . 4605 

1.4605 

1.4605 

1.4605 

.64 

1.5006 

1. 5006 

1.5006 

1.5006 

1.5006 

1.5006 

1.5005 

1. 5005 

1.5005 

1. 5005 

1. 5005 

. 66 

1.5406 

1.5406 

1.5406 

1.5406 

1.5406 

1.5406 

1.5406 

1. 5405 

1.5105 

1. 5105 

1. 5405 

.68 

1.5806 

1.5806 

1.5806 

1. 5806 

1.5806 

1.5806 

1.5806 

1.5805 

1.5805 

1.5805 

1.5805 

0. 70 

1.620 G 

1.6206 

1.6206 

1.6206 

1.6206 

1. 6206 

1.6206 

1. 6205 

1.6205 

1.6205 

1.0205 


TABLE I. COORDINATES OF Mx—x„, y — rjo ) — Concluded 
(b) Values of (i— *„)/S 


v-v, 
\ s 

'fx 

0 

0. 025 

0. 050 

0. 075 

0. 100 

0. 125 

0.150 

0.175 

0. 200 

0. 225 

0. 250 

— 0. 40 
-.38 
-.36 
-.34 
-.32 

0. 0258 
.0293 
.0332 
.0377 
.0428 

0. 0060 
.0151 
.0217 
.0281 
.0346 










— 0. 30 

0. 0485 

0. 0414 










-.28 

.0551 

.0489 

0. 0219 









-.26 

.0626 

.0572 

.0367 









-.24 

.0710 

.0663 

.0496 









-.22 

.0808 

.0766 

.0625 

0. 0256 








-0. 20 

0. 0918 

0. 0882 

0. 0761 

0.0500 








-.18 

. 1046 

.1013 

.0908 

.0700 

0. 0148 







— . 16 

. 1192 

. 1163 

. 1071 

.0897 

.0571 







-. 14 

. 1362 

. 1336 

.1254 

.1105 

.0854 

0. 0317 






-.12 

.1559 

.1535 

.1462 

.1330 

.1123 

.0782 






-0. 10 

0. 1791 

0. 1769 

0. 1702 

0. 1685 

0.1405 

0. 1137 

0. 0685 





-.08 

.2068 

.2047 

.1985 

.1877 

.1717 

.1490 

.1160 

0. 0572 




-.06 

.2408 

.2386 

.2326 

.2225 

.2076 

.1873 

.1598 

. 1204 

0. 0463 



-.04 

.2837 

.2816 

.2756 

.2654 

.2509 

.2317 

.2068 

.1743 

.1290 

0. 0408 


-.02 

.3437 

.3414 

.3344 

.3229 

.3072 

.2871 

.2624 

.2321 

. 1942 

. 1432 

0. 0487 
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Positive directions of axes and angles (forces and moments) are shown by arrows 


Axis 


Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

> 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 


X 

X 


L 

Y *Z 

Roll 


u 

P 

TifLtArftl 

Y 

Y 

Pitching 

M 

Z-—*X 

Pitch. 

e 

V 

« 


Z 

Z 

YAwing. . 

N 

X >Y 

Yaw 

'f' 

ID 

r 







• • 



Absolute coefficients of moment 

/Tf ^ 

^*“ 26 <s ^”~qcS 

(rolling) (pitchmg) 


O-^ 

^’^~gbS 

(yawmg) 


Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D 

P 

pID 

V' 

T 

Q 


Diameter 
Geotoetric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 

Thrust, absolute coefficient 



P 


0. 


n 

T 

n 

pn^D*‘ 


, Q 

$ 

pn^D* 



■P 

Power, absolute coefficient ^p— p^zj)6 

V pW 
Pr? 

Efficiency 

Revolutions per second, rps 
Effective helix angle=tan“’^ 2 ^^) 


5. NUMERICAL RELATIONS 


1 hp =76.04 kg-m/s=550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph=0.4470 mps 
I mps=2.2369 mph 


llb=0.4536 kg 
1 kg=2.2046 IV 
1 mi= 1,609.35 m=5,280 ft 
1 m=3.2808 ft 



